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I. INTRODUCTION 



The ongoing and planned high-precision measurements of the anisotropy in the cosmic 
microwave background radiation (CMBR) [1] may have serious impact on our views about 
the very early Universe. The large-angular-scale anisotropy in the CMBR is most likely 
caused by cosmological perturbations with wavelengths of the order and longer than the 
present day Hubble radius In- It is reasonable to expect that so long- wavelength perturba- 
tions are "primordial", survived from the epochs when the Universe was much younger. The 
wavelengths of the perturbations have enormously grown up since the time of generation 
but other physical characteristics of the perturbations can still carry imprints of their origin. 
This cannot be said with the same degree of certainty about the relatively short-wavelength 
cosmological perturbations (unless they are gravitational waves) which could have been dis- 
torted and contaminated in course of their life by many physical processes occurring in the 
Universe. 

It is remarkable that the origin of all three possible types of cosmological perturbations, 
that is the origin of density perturbations, rotational perturbations, and gravitational waves, 
may be of purely quantum-mechanical nature. Cosmological perturbations can be treated 
as excitations in gravitational field. In case of gravitational waves, they are just excitations 
in gravitational field itself. In case of density perturbations and rotational perturbations, 
they are excitations in gravitational field which accompany excitations in matter. In the 
very distant past, the density and rotational perturbations were excitations in the primeval 
medium that was filling the Universe at that time. The quantum-mechanical generation 
mechanism of cosmological perturbations relies only upon the existence of their zero-point 
quantum fluctuations and the nonvanishing parametric coupling of the perturbations to 
the variable gravitational field of the homogeneous isotropic Universe. The strong variable 
gravitational field of the very early Universe played the role of the pump field. It supplied 
energy to the zero-point quantum fluctuations and amplified them. More precisely, the 
initial vacuum quantum state of each mode of the perturbations has been transformed, as a 



2 



result of the quantum mechanical Schrodinger evolution, into a multiparticle state known as 
a squeezed vacuum quantum state. The generated perturbations have formed a collection of 
standing waves. The gravitational field of each of the three types of quantum mechanically 
generated perturbations can affect the propagating photons of the CMBR and produce 
anisotropy in CMBR. 

It was already emphasized [2] that there is a significant qualitative difference between 
gravitational waves on one side and density and rotational perturbations on the other side, 
with regard to possibility of their quantum mechanical generation. Gravitational waves 
oscillate in the absence of external gravitational fields, and their appropriate parametric 
coupling to the pump field follows directly from the Einstein equations. The parametric 
excitation vanishes only if the cosmological scale factor obeys the equation a" /a = 0, that 
is when there is no any pump field at all, a(rf) = const, or when the coupling a' is time 
independent. Up to this exception, one can say that the quantum mechanical generation 
of gravitational waves (relic gravitons) is unavoidable [3]. As for density and rotational 
perturbations, they are perturbations in matter being accompanied by perturbations of 
gravitational field. The ability to support oscillations of density and/or rotation and the form 
of their coupling to the pump field depend on a particular model of matter and its energy- 
momentum tensor. The very possibility of the quantum mechanical generation of these 
perturbations is model dependent. Recalling Einstein's definition of two pillars supporting 
general relativity, one can say that the quantum mechanically generated gravitational waves 
are associated with the pillar made of marble, while density and rotational perturbations 
are associated with the other one. 

A particular sort of matter that have received much attention in the recent cosmological 
literature, especially in the literature on inflation [4], is one or another version of a scalar 
field. Scalar fields is a nice theoretical model that has been used in physics in many different 
studies. Whether the global scalar fields do really exist in nature and, if so, whether they 
couple to gravity in the way we want, is presently unknown. However, we will follow the 
modern tradition in theoretical physics which states that everything that is not forbidden 

3 



is allowed. One can at least guarantee that a sort of inflationary expansion is a typical 
feature (attracting separatrix) in the space of homogeneous isotropic solutions to the Einstein 
equations with certain scalar fields [5]. Scalar fields cannot support rotational perturbations 
but they can support density perturbations. 

Specifically we will study a scalar field ip(t 3 ) with the energy-momentum tensor 
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where V((p) is an arbitrary scalar field potential, comma denotes a partial derivative. This 
model of the primeval cosmological medium satisfies both conditions for the quantum me- 
chanical generation of density perturbations be possible. First, the field can obviously 
support free oscillations in Minkowski space-time. Second, the explicit form of the energy- 
momentum tensor (1) reflects the appropriate (minimal, the same as for gravitational waves) 
coupling of the scalar field to gravity which was chosen by our will. So, on general grounds 
and by analogy with gravitational waves, one can expect that some amount of density per- 
turbations might have been generated by strong variable gravitational field of the early 
Universe. The problem is to quantify this expectation and to derive the observational pre- 
dictions, as reliable and detailed as possible, including the expected variations in the CMBR. 

Scalar fields and scalar field perturbations is a very popular subject in the framework 
of inflationary cosmologies. So popular, that many believe that the inflationary type of 
expansion is conditioned by the existence of scalar fields and that the very possibility to 
generate perturbations quantum mechanically relies on the existence of the De Sitter event 
horizon. This is not so. Inflation, if understood as a statement about the behavior of 
the time dependent cosmological scale factor, and not about creating and resolving the 
particle physics paradoxes, is a phenomenon more general than one particular realization of 
it with the help of a scalar field. [The attitude toward the architype inflationary solution - 
exponential expansion — has changed over the years. Astronomers of the older generation 
were embarrassed with the De Sitter solution but tried to apply it for the explanation of the 
galaxies' red shifts and statistics of quasars in the most recent Universe. Cosmologists of our 



time take the exponential expansion as something almost proven but apply it to the very- 
remote stages of evolution, somewhere near the Planck time.] And the quantum mechanical 
generation of perturbations is a phenomenon more general and universal than such concepts 
as global scalar fields, event horizons, and inflation. If it turns out that the inflationary 
hypothesis contradicts observations, the quantum mechanical generating mechanism will 
not die together with inflation. There is little doubt, for instance, that the search for relic 
gravitational waves will continue, with may be larger emphasis on relatively short waves 
rather than on long waves [6]. And a test of the quantum mechanical origin of cosmological 
perturbations will be a test of their origin, not a test of inflation specifically. 

The generation of density perturbations in inflationary models governed by the scalar 
field (1) was a subject of discussion in many research and review articles, and books. If one 
consults the most recent literature, one can find that the current situation is often summa- 
rized in the following, or similar, words: "Exponential inflation predicts a scale-invariant, 
Gaussian spectrum of scalar fluctuations and a smaller amount of tensor fluctuations .... 
Other inflationary models, for instance power-law inflation predict spectra slightly tilted 
away from scale invariance." (See, for example, [7] and references therein.) The expected 
amplitudes of density perturbations are usually quoted in the following, or equivalent, form 
(see, for instance, [8] and references therein): 

/6p\ hOT _ rnK 2 H 2 (<p) 

[p) x 87T3/2 mip) \ > 

where the quantities on the right- hand- side are supposed to be evaluated "when the scale A 
crossed the Hubble radius during inflation" . The denominator of this expression depends on 
the derivative of the Hubble parameter and goes to zero in the limit of exponential inflation. 
Apparently, this formula says that the predicted amplitudes of the scale-invariant spectrum 
are arbitrarily close to infinity, and the amplitudes of nearby spectra are "slightly tilted 
away" from infinity. According to this formula, the amplitudes of density perturbations are 
many orders of magnitude larger than the amplitudes of gravitational waves, if the expan- 
sion rate is sufficiently close to the exponential inflation. The belief that the amplitudes of 
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density perturbations are larger, or much larger, than the amplitudes of gravitational waves 
is considered to be a strong prediction of inflationary models based on the scalar field (1). 
For instance, the author of Ref. [8] concludes: "An observation violating this condition at 
any scale would immediately rule out the general class of models we are considering". The 
expected contribution of density perturbations and gravitational waves to the quadrupole 
anisotropy of CMBR was also under study. The authors of Ref. [9] (see also [10] and refer- 
ences therein) say that "The ratio of gravitational wave (T) to energy-density perturbations 
(S) contributions to the CMB quadrupole anisotropy is predicted to be T/S = 21(1 + 7)", 
where 7, in that paper, is the parameter in the equation of state for matter governing the 
inflationary expansion, p = jp. According to this formula, T/S vanishes in the limit of 
7 = — 1, that is in the limit of strictly exponential (De Sitter) inflation. Apparently, this 
formula for T/S is based on the authors' assumption that the effectiveness of generation of 
density perturbations is the higher the closer the expansion law to the exponential inflation, 
and goes to infinity in the limit of 7 = —1. Apparently, this is why T/S goes to zero in 
this limit. The statements about density perturbations are sometimes characterized as such 
that have been "widely studied and there is broad agreement regarding both methods and 
results ..." [11]. 

I suspect that the present paper will not belong to that category of studies that enjoyed 
the "broad agreement" ; my conclusions are considerably different from what was described 
in the preceeding paragraph. I will be arguing that there is no linear density perturbations 
at all at the purely exponential (De Sitter) inflationary stage for models governed by the 
scalar field (1). Density perturbations can only arise as a result of violation of the purely 
exponential expansion and transition to the radiation-dominated stage. Regardless of how 
close to zero was the derivative of the Hubble parameter "when the scale A crossed the 
Hubble radius during inflation", the today's amplitudes of density perturbations are finite. 
The amplitudes of gravitational waves are typically a little larger than the amplitudes of 
density perturbations, at least in the long wavelength limit where spectra are smooth and 
have the power-law behavior. Correspondingly, the contribution of density perturbations to 
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the quadrupole anisotropy is never much larger than the contribution of gravitational waves. 
In fact, it is somewhat smaller in the limit of long waves. 

As for the statistical properties of cosmological perturbations and, hence, the statistical 
properties of the CMBR fluctuations caused by them, it was already emphasized [12] that 
they are determined by the statistics of quantum states being generated, namely by the 
statistics of squeezed vacuum quantum states. 

Since the conclusions of this paper are in disagreement with other publications, we will 
present detailed derivations (which could have been omitted otherwise) in order to make it 
possible for the interested reader to compare the present calculations with those of other 
authors. The structure of the paper is the following. In Sec. II we present the general 
equations for density perturbations. The equations are applicable for matter with the 
energy-momentum tensor of arbitrary form. They use only the defining property of density 
perturbations, namely that the perturbations are based on the scalar functions of spatial 
coordinates. In Sec. Ill we apply these equations specifically to the initial stage 
(i stage) of cosmological evolution which is assumed to be governed by the scalar field with 
the energy-momentum tensor (1). No assumptions about a particular form of the scalar 
field potential V(ip) or a particular (for example, inflationary) type of expansion are being 
made a priori. The time-dependent coefficients of the differential equations for the pertur- 
bations are expressed in terms of the scale factor (and its derivatives) only. This reflects 
the underlying interaction of the perturbations with the variable gravitational pump field. 
The determination of all unknown functions describing density perturbations is reduced to 
solving a single differential equation which is very similar to the equation for gravitational 
waves. The behavior of solutions during a more or less gradual transition from the i stage 
to the radiation-dominated stage (e stage) is studied in Sec. IV. In order to deal with simple 
exact solutions at both stages we will be interested in a sharp transition from the i stage to 
the e stage. In Sec. V we apply the perturbation equations to the perfect fluid matter with 
arbitrary velocity of sound. We present solutions to these equations at the e stage and the 
matter-dominated stage (m stage) in the form convenient for matching the solutions at all 
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three stages (i, e, m). In Sec. VI we join the solutions, find the coefficients which were un- 
determined so far, and express the solution at the m stage entirely in terms of the functions 
(and their first time derivatives) describing the perturbations at the time of joining the i 
and e stages. As a preparation for quantization of the perturbations, we briefly discuss the 
density and rotational perturbations of matter placed in the Minkowski space-time, that is 
neglecting gravitational fields (Sec. VII). The quantization of density perturbations is per- 
formed in Sec. VIII. This procedure essentially repeats the steps which have been previously 
done for gravitational waves and rotational perturbations [12,2]. The quantum mechanically 
generated perturbations are placed in squeezed vacuum quantum states. Classically, one can 
think of the perturbations as of a stochastic collection of standing waves. The justification 
and necessity of the so-called Sakharov's oscillations in the spectra of density perturbations 
is discussed. In order to get the analytic results as detailed as possible, we specialize the scale 
factor of the i stage to the 77-time power laws which include inflationary models. This allows 
us to derive concrete power-law spectra of density perturbations at the m-stage. In Sec. IX 
we derive an exact formula for the angular correlation function of the CMBR temperature 
variations 5T/T caused by squeezed density perturbations. The multipole decomposition of 
the correlation function begins from the monopole term. The contributions to the monopole 
and dipole terms produced by individual waves with wavelengths exceeding Ijj are strongly 
suppressed, which is in agreement with previous results [13]. Nevertheless, one should be 
aware that not only the entire quadrupole but also some little portions of the measured 
mean temperature of CMBR and its dipole variation may be caused by perturbations of 
quantum mechanical origin. For one and the same cosmological model, the contribution of 
density perturbations to the quadrupole anisotropy of CMBR is never much larger than the 
contribution of gravitational waves, at least for models considered here. 
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II. GENERAL EQUATIONS FOR DENSITY PERTURBATIONS 



The unperturbed spatially flat FLRW (Friedmann-Lemaitre-Robertson- Walker) cosmo- 
logical models are described by the metric 

ds 2 = -c 2 dt 2 + a 2 (t)(dx l2 + dx 22 + dx 32 ) = -a 2 (r])(dr] 2 - dx l2 - dx 22 - dx 32 ) . (2) 

The scale factor a(rj) is governed by matter with the unperturbed values of energy density 
e , T° = -e , and pressure p , = poSf: 



3 (a'\ 2 



a 2 \ a 



K6 

2- 



2 " + " 

U V a , 



npo (3) 



where n = 87rG/c 4 and a prime is d/dr], d/dr] = (a/c)d/dt. The Hubble parameter is 
H = a/a = ca' /a 2 . 

It is convenient to introduce two new functions of the scale factor: 

a' a' 

<x(v) = ~ , l(v) = 1 - — • ( 4 ) 
a a z 

In terms of t-time the function 7 is 7(2) = —(H/H 2 ). Due to Eqs. (3) one has 

2a 2 

«(eo+Po) = — 7 • (5) 

(X 

The function 7(77) becomes a constant if a{rj) is governed by matter with the effective 
equation of state p = qe , where q = const. The scale factor a{rj) takes on the rj-time power 
law behavior 

a = l o \v\ 1+0 (6) 

(77 must be negative for expanding models with < 0) where l Q and (5 are constants. The 
constant l Q has the dimensionality of length. It follows from Eqs. (3) and (4) that 
7 = (2+/3)/(l+/3), q(/3) = (1 — /3)/3(l+/3), where the parameter (5 can vary in the interval 
— 00 < (3 < 00. In particular, 7 = 2 at the radiation-dominated stage, and 7 = 3/2 at 
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the matter-dominated stage. Note that 7 = in case of purely exponential (De Sitter) 
expansion for which a(t) ~ e Ht , H = const, a(rf) = lo^' 1 , j3 = —2. 
One can also derive from Eqs. (3) the relationship 

^ = -l + \l-^ = -U^aa^)'. (7) 
e 3 3«7 3a 

The ratio p' /e' becomes a constant for the scale factors (6), namely: p' /e' = q(f3). In 

particular, p' /e' goes to —1 in the limit of f3 — —2. 

The construction of density perturbations [14,15] (see also [16]) is based on the scalar 

functions Q^x 1 , x 2 , x 3 ) satisfying the equation 

Q^ i + n 2 Q = (8) 

valid in three-space dl 2 = dx l2 + dx 2,2 + dx s2 . For each wave vector n, one can choose two 
linearly independent solutions to Eq. (8) in the form e mx and e~ mx . From a given scalar 
field Q one can construct a vector field and two tensor fields: 5n~Q and Q^k = —niUkQ- 
For each n, the general perturbation of the energy-momentum tensor and the accompanying 
perturbation of the gravitational field can be written as a sum of products of time-dependent 
amplitudes and spatial functions introduced above. 

Without resticting in any way the physical content of the problem, it is convenient to 
work in the class of synchronous coordinate systems. (At this point the reader may have 
to be ready to exhibit certain resistance to the pressure from the proponents of the "gauge- 
invariant" formalisms.) Using the 77-time coordinate, one can write the general expression 
for the metric tensor including perturbations as: 



ds = —a 1 



drj 2 — (Sij + hij)dx l dx^ 



g 00 = -a 2 , 9oi = 0, gij = a 2 (1 + hQ)Sij + hn 2 Q )id . (9) 

The function h{rj) represents the scalar (proportional to SyQ) perturbation of the gravita- 
tional field while the function hi(rj) represents the longitudinal-longitudinal (proportional to 
riirijQ) perturbation. The general expression for including perturbations can be written 

as 

10 



a 2 a 2 cr 

T t k = Vo 5l + l(p! + Pl )Qtf + ^n~ 2 W Q/ . (10) 

The form of Eqs. (9) and (10) is based solely on the definition of density perturbations 
and our choice of synchronous coordinate systems. In all other respects, the representation 
(9) and (10) is general. The particular notations for arbitrary functions describing the 
perturbations are chosen for later convenience. 

The arbitrary functions h(rj), hi(rj), ei(rj), pi(rj) pi(rj), £'(r)) should satisfy all together 
the perturbed Einstein equations: 

3ah' + n 2 h — ah[— K6\ (11) 

h'= <' (12) 

-h" - 2ah'= k Pi (13) 

^(h' l ' + 2ah' l -n 2 h)= m . (14) 

There are too many unknown functions to be found from Eqs. (11)-(14). This requires 
us to specify a model for matter and its energy-momentum tensor. We will consider three 
consecutive stages of expansion: i stage governed by the scalar field (1), and the subsequent 
e and m stages governed by perfect fluid with the energy-momentum tensor 

T/ = (e + p)«X+PV- ( 15 ) 
The equation of state at the e and m stages is p = |e and p — 0, respectively. 
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III. DENSITY PERTURBATIONS AT THE INITIAL STAGE OF EXPANSION 

GOVERNED BY A SCALAR FIELD 



At the i stage, the evolution of the scale factor 0(77) is determined by the unperturbed 
homogeneous scalar field ip = <p (jf)- The unperturbed values eo, po are given by Eq. (1): 

eo = ^' ) 2 + V( V ) (16) 
Po=^W-V(<p) ■ (17) 

By summing up Eqs. (16) and (17) and comparing the result with Eq. (5) one can derive 
the equation 

«K) 2 = 2a 2 7 . (18) 

It follows from this equation that 7 > for the scale factors governed by the scalar field (1). 
The De Sitter case corresponds to <p' = 0, ip = const and e = — po = V(<po) = const. 
If ip' Q 7^ 0, one can use the equation 

e' = -3a(e + p ) 

which is a consequence of Eq. (3), and obtain with the help of Eqs. (16) and (17): 

<pl + 2a<p' + a 2 V^ = (19) 

where V jtp = dV((p)/dip, the derivative is taken at ip — ip . The further useful relationships 
following from Eqs. (18) and (19) are: 

ip a 2 7 

-—V >ip =2a+- + -+- . 21 
^0 a 2 7 

The perturbations of the gravitational field are associated with the perturbations of the 
scalar field. We will write the perturbations of the scalar field as 

up = Mv) + <Pi(v)Q ■ (22) 
12 



Having at our disposition the energy-momentum tensor (1) and the definitions (22), (9), 
(10) we can directly calculate the functions e±, pi, pf. 

ei= </?o</?i + a?(piV v (23) 

e= -y^i (24) 

Pi= <fo<p'i ~ a ViKv ( 25 ) 
Pi=0 ■ (26) 

We will now assume that <p' ^ 0. The De Sitter case y?g = will be considered separately 
at the end of this Section. It follows from Eq. (24) that (p 1 = — £'/Vo- Inserting this value 
of ipi into Eqs. (23) and (25) and using Eqs. (20) and (21), one can express ei, pi, in terms 
of h(rj): 

K€l = -h" + ti (2a + 2^ + ^ (27) 
KPl = -h" - 2ati . (28) 

We should now return to the perturbed Einstein equations (11)-(14) making use of 
Eqs. (27) and (28). Equation (11) can be written as an expression for h'^rf) in terms of 

1 



K 



h " + h i[ a -2--^)+n 2 h 



a \ a 7 



(29) 



Equation (12) expresses £' in terms of h! . Equation (13) is satisfied identically. Equation (14) 
reads as 

h" + 2ah' l -n 2 h = . (30) 

Thus, if one knows the function h(r}), all other functions describing the density perturbations, 
namely hi{rj), Pi(r)), ei(r)), and £'(77), can be found with the help of Eqs. (29), (28), (27) 
(or, equivalently, (11)) and (12). To derive the equation for h{rj) we substitute Eq. (29) into 
Eq. (30) and obtain 
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ti" + h"Ua 1 -'£\ +h' 



n*-2a' + 2 1 a*-- 1 --r- 

a 7 V 7 



+ n 2 a-fh = 0. (31) 
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Equation (31) is a third-order differential equation. There should be no wonder (and no 
panic) on this occasion. One of solutions to this equation we know in advance, this solution 
is 

h = C- (32) 
a 

where C is an arbitrary constant. We could have expected the existence of this solution, even 
before solving the equation for h(i]) , because the perturbation of this form can be generated 
by a coordinate transformation which does not violate our choice of synchronous coordinate 
systems and, hence, does not destroy our initial form (9) of the perturbed metric. (One can 
easily check that the function (32) is indeed a solution to Eq. (31).) Concretely, one can 
perform a small coordinate transformation 

In terms of new coordinates fj, x l the transformed components (9) take on the form 

goo = -a 2 (fj) , g i = , g ik = a 2 (f]) (1 + hQ)5 ik - hin~ 2 Q jijk 
where 

h = h + C-, h = h { + Cn 2 I a' 1 d V . (33) 
a J 

The same transformation should be applied to the components of the energy-momentum 
tensor. Even if the original h, hi are zero, the transformed h, hi are not zero. After erasing 
the overbars in the transformed components of the metric, one returns to Eq. (9). The 
freedom of choosing different freely falling coordinate systems and corresponding spatial 
slices 7] = const gets represented in the form of freedom to choose different solutions from 
the family of all solutions for the perturbations. All choices of C are equally well "physical" . 
The integral in Eq. (33) produces an additional integration constant which reflects the 
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possibility to make a purely spatial transformation and to shift hi by a constant value, but 
we will not actually need this coordinate freedom. It follows from Eq. (33) that there are 
two functions (and many algebraic and differential combinations constructed from them) 
that do not contain C at all: 

u = h' + ajh , v — h', n 2 h . 

a 

The solution (32) allows one to reduce the third-order differential equation (31) to the 
second-order differential equation. To reach this goal we use the function u(rj): 

u = ti + ajh . (34) 

Obviously, this function vanishes on the solution (32). By substituting Eq. (34) into Eq. (31) 
we derive the equation for u{rj): 



u" + u [ 2cry — — ] + u 
V 1 



n —2a — a I — I 



. (35) 
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Note that the coefficients of this differential equation depend exclusively on the scale factor 
and its derivatives. This fact is a manifestation of the underlying interaction of the pertur- 
bations with the cosmological pump field. No special assumptions about the shape of the 
potential V(ip) or such things as "nonsimultaneous rolling the scalar field down the hill" 
have been made whatsoever. 

Our next move is to transform Eq. (35) to the form similar to the equation for gravita- 
tional waves. This will allow us to use certain results derived previously for gravitational 
waves and rotational perturbations [17,2]. In order to get rid of u' we introduce the function 
/i(r)) according to 

u = — - — n . (36) 

CL 

It follows from Eq. (18) that the function 7(77) is nonnegative if the scale factor is governed 
by the scalar field (1) which we study here. However, Eq. (35) is formally applicable to 
negative 7 as well. It may happen (as the author thinks) that Eq. (35) has a wider domain 

15 



of validity and can be used, for other models of matter, with negative 7 too. If this is the 
case, one is free to modify Eq. (36) by using instead of ^7. Anyway, with the help of 
Eq. (36) one derives the equation 

H" + fi[n 2 - U(rj)} = 



where 



U(r]) = a 2 + a +a— + 



/, 7 , 1 ( 7 \ , 1 { 7 \ 



+ 



7 4y7y 2 y 7 
a" „ , , 1 



U ( V ) + UM 



U (r})= a 2 + a' = — , ^(77) = — 



a 



ail > - I 7 ' 2 + I 77 " 



(37) 



The effective potential U{rj) can also be written as U(rj) = {a^Jl)" / a^Jl which reduces our 
basic equation to the form 

.2 (av/77" 



n" + n 



n 



= . 



(38) 



(The function a^/7 can be related with the function z discussed in [18], see also the early 
papers [27].) 

Let us recall [3] that in the case of gravitational waves the potential U (rj) consists only 
of the Uq(t]) term, so that the basic equation is 



n 2 - 







(39) 



The gravitational wave potential U depends only on the first and second time derivatives 
of the logarithm of the scale factor: (In a)', (In a)". The potential U(rf) for density per- 
turbations is more complicated and includes also (In H)', (In H)", and (In H)'". We note, 
however, that the potentials are exactly the same, and, therefore, the basic equations and 
solutions for density perturbations and gravitational waves are exactly the same, if 7 is 
constant, that is for the scale factors (6). For this class of pump fields, the general solution 
to Eq. (39) can be written in the form (for non-half-integer (3): 



fi(ri) = {nr]) 1/2 



AiJ^Unrj) + A 2 J_ {p+ i } (nrj) 



(40) 
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Having two linearly independent solutions to Eq. (38) one can construct h(r)) and, hence, 
to find the rest of functions describing density perturbations. It follows from Eqs. (34) and 
(36) that 

U = - ia h + (41) 
a 

and 

Kv) = - [\Vl dv + -Q (42) 
a Jvq a 



where r] Q is some initial time where the initial conditions are to be imposed. The constant C 
entering Eq. (32) is denoted Cj at the i stage and will have the labels e and m at the e and 
m stages. All (complex) solutions to our perturbation problem for a given wave vector n 
are completely determined by three arbitrary and independent (complex) constants. Two of 
them define a solution to Eq. (38) (these constants are A 1 , A 2 when Eq. (40) is applicable). 
The third constant, Cj, describes the remaining freedom in our choice of coordinates. This 
remaining freedom is not a misfortune of the theory. On the contrary, it will later allow us 
to join our coordinate system right to the comoving synchronous coordinate system at the 
m stage. 

One can show by using Eqs. (24), (12), and (41) (and assuming ip' ^ 0) that 

1 



^i(^) = "i= -v- \h h 

V 2k [a 

One can also find with the help of Eqs. (11), (13), and (41) that 



(43) 



Pi 



^ / A ■ (44) 

H'-nUa+^\ +3a^fah 

The quantity q, where cf/c 2 = pi/ei, plays the role of the velocity of sound for the high- 
frequency scalar field oscillations (see also Sec. VII). 

Similarly to what is true for gravitational waves, solutions to Eq. (38) are different 

for the high-frequency and low-frequency regimes. In the former case, n 2 ^> |f/(^)| and 
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pi e ±mv . In the later case, n 2 <C \U(rj)\ and two independent solutions are /ii ~ 0^/7, 
ji 2 ~ a ^/l 5 dv/i a ^) 2 - The functions fj,i, 112 generalize the corresponding solutions for 
gravitational waves by replacing a with Specifically for the scale factors (6), the 

solutions /ii, jj,2 are /ii ~ i] 1+l3 and fi 2 ~ in agreement with Eq. (40). 

In the high-frequency regime, the term p! dominates the other terms in Eq. (44). As one 
could expect, in this regime, the velocity of sound is equal to the velocity of light, cf = c 2 . 
In the low-frequency regime, that is when a given mode enters the under-barrier region, the 
dominant solution is p\ (for a review, see [19]). Using this solution in Eq. (44) one can show 
that, in this regime, p\/e\ pa g(/3), that is the "velocity of sound" is the same as the one 
defined by p' /e' Q . In particular, pi/ei goes to —1 for the low-frequency scalar field solutions 
at the De Sitter stage, (3 = —2. 

Equations (31), (35), and (38) have been derived under the condition ip' Q 7^ 0. However, 
the final formula (42) gives the correct result h(rj) = Cj/Zo = const in the De Sitter limit 
7 = 0, (p' = 0. One can analyze this case separately, referring to the starting Eqs. (23)-(26). 
One can see that Eqs. (24) and (12) give £' = 0, h — const. Equations (13), (25), and (23) 
give pi — 0, V :lf> — 0, ei = 0. Finally, Eq. (11) requires h[ = —r\n 2 h = —r]n 2 Ci/l . But 
this solution for h, h[ is precisely solution (33) which can be eliminated by a coordinate 
transformation. In the De Sitter case governed by the scalar field (1) there is no density 
perturbations at all. Note that the function (pi(rj), Eq. (22), remains arbitrary and the 
wavelengths of these fluctuations are growing in the course of expansion. If one wishes, one 
can attach to tpi such words as "inflation is pushing the waves beyond the De Sitter horizon" . 
Nevertheless, the result will be zero, as long as ipi is not accompanied by perturbations of 
the gravitational field. This is an instructive example in order to realize that to "stretch the 
waves outside the causal horizon" is not all we need for generation of density perturbations 
(likewise, it is not sufficient to simply stretch the electromagnetic waves "beyond the horizon" 
in order to generate photons). 
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IV. LATE TIME EVOLUTION OF THE PERTURBATIONS AT THE INITIAL 

STAGE 



The main uncertaintes about the evolution of the very early Universe refer to the times 
that preceded the epoch of the primordial nucleosynthesis. Whatever was the initial stage, it 
supposedly went over by that epoch into the radiation dominated stage governed by the scale 
factor a(r)) = l a e (v — Ve)- The constants a e , r\ e are to be determined from the continuous 
joining of a(rj) and a'(rj) at the time 77 = 771 of transition from the % stage to the e stage. If 
the % stage is described by the scale factors (6), one derives 

In further applications, we intend to use simple solutions (40), (42) and to make their 
appropriate joining with perturbations at the e stage. However, the function 7(77), being 
equal to the constant 7 = (2+/3)/(l+/3) at the i stage, and to the constant 7 = 2 at the e 
stage, experiences a finite jump at the transition point 77 = 771 . This presented no problem for 
gravitational waves, since 7' (77) did not enter the gravitational wave potential Uo(i]). But this 
becomes important for density perturbations, since the Ui(t]) part of the potential acquires 
increasingly growing values at the end of the i stage for steeper and steeper transitions. 

To deal with the problem, we introduce a parameterized set of smooth functions 7(77) 
that approximate the step function in the limit of the parameter e going to infinity: 

, \ 4 + 3/3 , r . . A . 

^ ) = WW) + WW) & ■ ( 5) 

For large negative values of 77 the function (45) goes to (2 + /3)/(l + /3), and for large positive 
values of 77 it goes to 2. We may surround the transition time 77 = 771 by a thin "sandwich" 
with boundaries at 771 — a and 771 + a. The asymptotic values of 7(77) are already reached 
with arbitrary accuracy at the boundaries, if e is sufficiently large, e 3> 1/er. 

The function (45) can be integrated, see Eq. (4), to produce the function a(rj): 
1 77 



a(ri) 1 + p 2(1 + p) 



I e t(v-vi) _|_ e -t(v-vi) 
77 - 7/! + - In 
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(46) 



Again, for sufficiently large e, the function a(r)) quickly approximates (1 + f3)/r) to the left 
of the transition point, and 1 / (77 — i] e ) to the right of the transition point. These are the 
values of 01(77) that are appropriate for the % stage (6) and the e stage, respectively. 

The divergent functions 7', r )' 2 and 7" participate in the potential Ui(rf) (Eq. (37)). The 
function 7' grows as e at the point r\ — r\\. The function 7" is equal to zero at 77 = 771, but 
it grows as e 2 slightly to the left of this point, and it grows as — e 2 slightly to the right of 
this point. We assume that the transition to the e stage has completed at 77 = 771 + a, and 
we let a to go to zero. The function Ui(rj) compresses and stretches to the arbitrarily large 
positive and negative values when a goes to zero and e goes to infinity. Examining Eq. (38), 
one can expect that the value of 7/(77) at 77 = 77! + a will be different from the value of 7/(77) 
at 77 = 771 — a. The integration of //' in the limits from 771 — a to 771 + a gives a jump in 7/ 
which depends on the value of the integral from the divergent part of the potential Ui(rj). 
Fortunately, it is not 7/(77) itself, but a particular combination (^/a)[fi' — fi(a + 7Y27)] 
that we will need to know in our further calculations. This simplifies the analysis. Due to 
Eqs. (29) and (41) this combination is precisely the function ^(77) introduced in Sec. III. 

In terms of the function v (77), where 

?j = — fi -fi[a+-—\ 
a I \ 2 7 / 

the basic equation (38) takes on the form 

(a 2 v)' = -n 2 a^fi . (48) 

The integration of this equation over the thin "sandwich" shows that v\ v=Vl+0 = v\ v=m - . 
In other words, the function fi' — 71(0 + |^r) taken right at the beginning of the e stage is 
equal to the value of this function taken right at the end of the % stage (6) times the factor 

7f Vl+f • In addition to the conditions: 7^=^-0 = y^fjf , 7^=^+0 = 2, /jl\ v = Vi -o = v\ v = m +o, 
this establishes the rules for going through the "sandwich" . 
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(47) 
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V. DENSITY PERTURBATIONS IN THE PERFECT FLUID MATTER 



We will now consider Eqs. (11)-(14) at the perfect fluid stages governed by the energy- 
momentum tensor (15). Similarly to the scalar field case, the longitudinal-longitudinal part 
of stresses vanishes, pi = 0. For the easier handling of arbitrary ei, p\ one can introduce 
the following notations: pi/ei = c 2 /c 2 and p' /e' = c 2 s /c 2 , see Eq. (7). These definitions are 
convenient but, generally speaking, they have only formal meaning, since both pi/ei and 

Pi/ 

e' can be negative. However, in certain regimes, the quantity q is a genuine longitudinal 
velocity of sound (see Sec. VII). Our first intention is to derive the equation for h(r/), anal- 
ogous to Eq. (31) and valid for arbitrary nonzero p\je\. Specific and cf = 
will be considered separately. 

In order to derive the equation for h(i]) one can essentially repeat the steps that have 
lead to Eq. (31). Find h'^rf) from Eq. (11) and plug it into Eq. (14). Use Eq. (13), the 
first derivative of this equation, and the definition of cf/c 2 . As a result, one arrives at the 
equation 



h"'+h" 
+ ti 



-2y 



3a + «7 + 3a— + 



2H 



(<?)' 



c 

+ hn 2 a^-^ = . 



(49) 



Now, introduce the function u(rj) according to Eq. (34) and use it in Eq. (49). Equa- 
tion (49) can be reduced to 



u"+u' 



+u 



3a + 3a% - 

c 2 cf 

n 4 +(<*+*) Li 



a 



c 2 

2_, s 



+hl 3a 7^ 2 



3a 



„2y 



a" a" 
+ — + 2— 



a 



+ 



a 

= 0. 



(50) 



The last term in this equation vanishes if 

(cf)' 



c 2 
3a — 



r 2 



C 2 

3a% 



(c 2 s y 



(51) 
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which integrates to 



c 2 s const 
cf aa 2 j 



An assumption which is usually made for perfect fluids is: 

cf = c\ . (52) 

Note that Eq. (52) is certainly true for matter with the equation of state p = qe, where q is 
a constant, but Eq. (52) is not true in general, and it is not true for the scalar field matter 
(1) (unless one considers the under-barrier region where Eq. (52) is true approximately, see 
Sec. III). Due to Eq. (51) the last term in Eq. (50) cancels out. (If we have not assumed 
(51), the function h{rj) = C(a/a) would not have been a solution to Eq. (49).) 
We can now introduce the function u{rj) according to (compare with Eq. (36)): 

U = — - — C s V . (53) 

a 



In terms of v(rj), Eq. (50) takes on the form 



v + v 
where 



_,2 



n 2 % - W( V ) 



(54) 



a a a^^Cs 

The potential W(rj) depends exclusively on the scale factor and its derivatives. Having a 
solution to this equation and using Eqs. (53) and (34), one can construct the function 
h(rf) and the rest of perturbations. Similarly to the scalar field case, all solutions for pertur- 
bations with a given n are defined by three constants one of which describes the remaining 
coordinate freedom. These constants are expressible in terms of the constants given at the 
preceding i stage through the joining of the perturbations at the transition time i] = rji from 
the % stage to the perfect fluid stage. 

We will now consider Eq. (54) specifically at the radiation-dominated stage p = |e. We 
have cf = c\ = (l/3)c 2 , 7 = 2, and the scale factor 
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a(rj) = l a e (r] - i] e ) (55) 

where the constants a e , r\ e are to be determined from the continuous joining of a{rf) and 
a' (77) at the transition time 77 = 771 . The potential W(rf) vanishes. Equation (54) simplifies 
to the familiar equation 

v" + \n 2 v = (56) 

which describes the time- dependent part of sound wave oscillations in the radiation- 
dominated fluid. In what follows, we will be using the general solution to this equation 
written in the form 

v = B l e -^ {T '-" e) + B 2 e*^-*) (57) 

where Bi, B 2 are arbitrary and independent (complex) numbers for each individual wave 
vector n. 

The function h{rj) is determined by a known solution for 1/(77) an d a coordinate solution 
with arbitrary constant C e : 

h( v ) = - r udrj + -C e . (58) 
a Jrii a 

All other functions are expressible in terms of h{rj). In particular, 

h' = -[3h" + 9ah' + n 2 h] . (59) 
a 

The general solution (57) is always oscillatory in 77-time. The perturbed energy density, 
pressure, and the associated gravitational field h(rj), hiirf) oscillate in space and time as 
they should do for sound waves. However, if one considers these oscillations at intervals of 
time shorter than their period, they may appear as consisting of "growing" and "decaying" 
solutions. In particular, this happens if one considers relatively long waves, 
75(^1 ~ Ve) = 2/i <^ 1, at their first oscillation since the beginning of the e-stage, that is 
while the condition -^(v — i] e ) <^ 1 is satisfied. Under this condition, the function h(r)), 
Eq. (58), can be approximated as 
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where 



— (V-Ve), E 1 - Bl + B \ B, ~ m{Bl - B2) 



and 



ol a e [ n 

The common practice [15] is to use the coordinate freedom for elimination of the "most 
divergent" term in Eq. (60), which is also the "most decaying" term, if one goes forward in 
time. This is achieved by such a choice of C e that C e = 0, and the first term in Eq. (60) 
vanishes. Then, the energy density perturbation 5e/eo (use the definition 

Se/e = —Q 

and calculate ke\ according to Eq. (11)) can be approximated as 

| = ^B.n + l -n 2 B 2 rf + ■ ■ -)q . (61) 

The part of Eqs. (60) and (61) which depends on the coefficient B\ is usually called the 
"decaying" solution, while the part with the coefficient B 2 is called the "growing" solution. 

Despite the possibility of identifying (quite artificially) the "growing" and "decaying" 
solutions, density perturbations at the e-stage form a collection of traveling sound waves 
with arbitrary amplitudes and arbitrary phases, as long as constants B±, B 2 are arbitrary 
and independent. One should not think that simply because the sound waves have spent 
some time "beyond the horizon", nfj 1, they would transform into standing waves at later 
times of their history after they came "inside the horizon", nfj 1. To illustrate this point, 
let us take into account the spatial part of the perturbations and consider the contribution 
/i n (?7,x) of a given mode n to the total field h(r],x.) = S n ^n(^, x )- This contribution can 
be written as 
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h n (r], x) ~ hn e mx + h* n e" inx ~ (fi ln e~ infi - B 2n e inf >)e in * + (#*„ e m ^ - 5 2 * n e -^) e - inx 

= 2\B ln \ cos(nfj - nx - (p ln ) - 2\B 2n \ cos(m? + nx + (p 2n ) (62) 

where Bi n = iBi^e^ 111 , B 2n = \B 2n \e l(fi2n . The last line in Eq. (62) shows explicitely that, in 
general, one is dealing with waves traveling in opposite directions with arbitrary amplitudes 
and arbitrary phases. A standing wave can only be "generated by hand" , by assuming that 
the constants B in , B 2n are strictly related. This happens if one declares that he/she is only 
interested in the "growing" solution and puts B ± — 0. Then, the complex amplitudes B lri , 
B 2n become related: |B 2n | = (— l) fe+1 |i?i n |, (p 2n = (p ln — kir, (k = 0, 1,2, • • •), and the last 
line in Eq. (62) can be transformed to 

h n (r], x) ~ 4 1 .Bi n | cos nfj cos(nx + (pi n ) 

which is a standing wave indeed. 

Standing sound waves at the e stage are responsible for so called Sakharov's oscilla- 
tions [20] in the power spectrum of density perturbations in the present Universe, at the m 
stage. As we have shown, standing sound waves cannot originate somehow automatically at 
the e stage, simply because of the transition from the "growing" / "decaying" regime to the 
oscillating regime (see also [21]). If one works with classical density perturbations at the 
e stage and makes no additional assumptions, one can say nothing about the necessity of 
standing waves, except of postulating this. The point of this discussion is that the quantum- 
mechanical generating mechanism, which we are considering in this paper, does really create 
standing waves. Standing waves arise for gravitational waves, rotational perturbations, and 
density perturbations. The physical reason for this is that the waves (particles) are gener- 
ated in correlated pairs with equal and oppositely directed momenta (the two-mode squeezed 
vacuum quantum states). This is true for waves of any wavelength, as soon as conditions for 
their generation are satisfied. Technically, as we will see later, the second term in Eq. (60) 
taken at the beginning of the e stage turns out to be much smaller than the third term in 
Eq. (60), that is B x + B 2 w for Vl < 1. 
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We should now discuss a great difference between sound waves and gravitational waves 
with regard to their evolution in time. The velocity of sound waves at the e stage is only 
times smaller than the velocity of gravitational waves, but their amplitudes behave 
drastically different. The amplitude of a gravitational wave decays as a" 1 in course of time. 
The amplitude of a sound wave, as one can see from Eq. (58), decays as a~ 2 , since a ~ a' 1 . 
This leads to a difference in solutions even for relatively long waves which did not complete 
even one cycle of oscillations during the entire e stage from rj = rji to r\ = r/ 2 - As an 
illustration, let us consider sound waves which barely reached the oscillating regime by the 
end of the e stage. Their wave numbers satisfy the condition -^(i]2 — i] e ) = ^- = y 2 ~ 1. 
These are the waves whose wavelength was of the order of the Hubble radius at the time of 
transition from the e stage to the matter dominated m stage. Assuming that the present 
day Hubble radius is Z# ~ 6 • 10 3 Mpc and that the present day scale factor a(r] R ) is 
a(rj R ) ks 10 4 0(772), their wavelength today A c = 2ita(riR)/n c is about 220 Mpc. The usual 
practice, in addition to eliminating C e , is to concentrate on the "growing" solution, that is 
to assume that at the beginning of the e stage the second term in Eq. (60) is smaller, or 
at least not larger, than the third term. Under these conditions, the final numerical value 
of the function h{rj) is of the same order of magnitude as the initial value, h{rj 2 ) ~ h(r)i), 
for the wavelengths of our interest. In other words, if the preceding i stage produced h(r)) 
with some initial numerical value h(rji), this number will effectively be transmitted to the 
beginning of the m stage. 

However, in the very same coordinate system where the "most decaying" term in Eq. (60) 
was eliminated, the time derivative h'(rj) was left large. The final value of h'{rj) for the waves 
of our interest is ^'(772) ~ nhfa). According to Eq. (12), the function h'(rj) describes velocity 
of matter. This velocity will be inheritted by matter at the matter-dominated stage. But this 
is not velocity of the fluid elements with respect to each other, this is not velocity describing 
deformations of the medium, and the functions h'(rj), h'^rj) are not the ones that we may 
use for our later calculation of the varations in CMBR. The function h'(rj) describes velocity 
of matter with respect to the coordinate system that we have chosen for our convenience of 
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eliminating the "most decaying" term. At the m stage, however, we are more interested in 
the comoving coordinate system. We are interested in the density contrasts and deformations 
of the fluid itself, we are interested in the components of the accompanying gravitational 
field that we may use for calculations of 5T/T. We need a coordinate system which provides 
vanishing of h'(rj) by the end of the e stage. This requires a different choice of C e . Under 
this new choice of C e , the first term in Eq. (60) survives, and numerically the same, as in 
the previous example, initial value h(r)i) transforms into a small number ^(772) ~ Z/iM^i) by 
the beginning of the m stage. This is what we need to keep in mind when we will compare 
the amplitudes of density perturbations and gravitational waves. 

Finally, let us consider density perturbations at the m stage, p — 0. At this stage, one 
has cf — c 2 s — 0, 7 = §, and the scale factor 

a(rj) = l a m (rj - r] m ) 2 . (63) 

The scale factor and its first time-derivative are continuous at the time r\ = 7? 2 of transition 
from the e stage to the m stage. Therefore, a rn = a e /4(r? 2 — r) e ), i] m = — r? 2 + 2i] e . It follows 
from Eqs. (13) and (14) that the general solution for h(r)), hi(r)) has the following form 

h=C\-\ C m , h — —— C m , 
a la 

h[ = \c x n\r, - Vm ) + -n 2 C m + C 2 { f ~ ^ . (64) 
o a [rj — r] m ) 

The energy-density and velocity perturbations can be found from Eqs. (11) and (12). Simi- 
larly to the preceding % and e stages, the perturbations are completely determined by three 
constants C±, C2, C m one of which, C m , reflects the remaining coordinate freedom. 

The constant C m is entirely responsible for a possible relative velocity of our fluid with 
respect to a chosen synchronous coordinate system, see Eqs. (12) and (10). A given coordi- 
nate system is not comoving, T l 7^ 0, as long as C m 7^ 0, h! ^ 0. But we know that for a 
dust-like fluid without rotation one can always introduce a coordinate system which is both 
synchronous and comoving. This is reflected in our ability to remove the C m term from h(rj) 
and hi{rj) by a coordinate transformation (33). Thus, the choice C m = in Eq. (64) is not 
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a restriction of the physical content of the problem, it is an allowed choice of the coordinate 
system. 

It is here that we will eventually restrict our coordinate freedom, we will put 

C m = . (65) 

The constants C e and C, will not be arbitrary any longer, they will be determined from 
the continuos joining of solutions. We need the comoving coordinate system for simple 
and appropriate formulation of the 5T/T problem. We are interested in the temperature 
of CMBR and its anisotropy seen by a comoving observer, that is by an observer whose 
world line is one of the matter's world lines. One of these idealized comoving observers is an 
observer on Earth (up to accuracy of some nonzero peculiar velocity, local rotation, etc.). 
We are much less interested in feelings of an observer who wanders in the Universe with 
arbitrary time-dependent velocity. In the comoving coordinate system, the world line of a 
comoving observer is described by simple equations x % = const. 
Upon the choice of C m = 0, the perturbations reduce to 

h(v) = Cx 

h^ri) = —Cm^rj - Vm ) 2 - -C 2 ~ Vm )l . (66) 
tyn 10 W ' ; 3 (r)-Vm) s 

From Eqs. (66) and (11) one can derive the familiar expression [15] 



6e 
eo 



lm) 3 



Q. (67) 



(One may wish to correct a misprint in Eq. (115.21) of Ref. [15]: the decaying solution 
behaves as r]~ 3 , not as the printed r]~ 2 .) As long as the constants Ci, C 2 are arbitrary, 
the power spectrum of the density perturbations is arbitrary. In particular, there is no 
Sakharov's oscillations, a priori, and they do not arise simply because of the transition from 
the e stage to the m stage. For instance, one can start from a perfectly smooth spectrum 
at i] = r/2 and extrapolate these data back in time up to the beginning of the e stage. 

For the further calculations of 5T/T we will need the first time derivatives of the gravi- 
tational field perturbations at the m stage in the comoving coordinates. Since for the scalar 
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component of the perturbations one has h! = 0, it is only the longitudinal-longitudinal 
component h[ that is effective. 
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VI. JOINING THE PERTURBATIONS AT THE THREE STAGES 



We are now in the position to start our operation of joining the solutions at i, e, and 
m stages. We want to derive from the first principles the expected density perturbations at 
the m stage. Of course, the result will depend on the unknown behaviour of a{rj) at the i 
stage. But this is precisely why we are doing this study: we try to learn something about 
the evolution of the very early Universe by deriving the expected variations in CMBR and 
comparing them with the observations. 

The general rule for joining solutions to Einstein's equations is to match from the both 
sides the intrinsic and extrinsic curvatures of the transition hypersurface [22]. For our 
solutions written in the class of synchronous coordinate systems, this translates into the 
continuity of the spatial metric and its first time derivative. Since we have already assumed 
that a(r]) and a'(rj) join continuously, it is the continuous joining of h(rj), hi(r)), h'(r)), and 
h'^rj) that should be ensured. In fact, it is sufficient to follow h(r)), h'ijf), and h'^rf) as hi(rj) is 
derivable from h'^rj) at all three stages up to the integration constant which can be removed 
by the remaining integration constant in Eq. (33) anyway. It is convenient to write h[ at the 
i and e stages, respectively, in the form (use Eqs. (29), (59) and (42), (58)): 



We will denote a(rj), a(r)) at rj = rji and r\ = 772 by ai, a± and 02, «2 respectively. It is also 
convenient to introduce the parameter y — -^{t} — r) e ) and its values y±, y2 at the transition 
points. 

We will first make the joining of solutions in general form, without adopting any particu- 
lar coordinate system, any particular behavior at the i stage, and any particular wavelength 
of the perturbations. 

Let us start from the i-e transition. Whatever was the % stage and its late time behavior, 
it produced certain ii(r)i), //(r/i) and ended at r\ = rji with some values of the scale factor and 




h\= —h + ~{v' -av) 



(68) 



a a 
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its derivatives. For generality, we do not assume, for the time being, that 7(771) is exactly 2 
and 7' (771) is exactly zero. 

From the continuous joining of h(rj), h'(r)), and h'^rj) one derives: 



J no 



B 1 e- tyi + B 2 e iyi = Vt> + a(2 - 7) 



rvi 
/ 



+ c 



5! e^ 1 (1 + i Vl ) + 5 2 (1 - i Vl ) = --aV7 



/i - jjL [a + 



If 

27, 



(69) 



(70) 



(71) 



All functions in these equations are taken at 77 = 771 so that, for instance, 7 means 7(771), fjf 
means //(771), 7' means 7 / ( r 7i) ) etc. For a more compact record we will also use the following 
notations: 



"l = — yi(m)v(vi) , v i = 

ttl CL\ 



I 1 7' 



and 



C — Cj + 



Jrio 



Equations (69)-(71) allow us to express the constants £>i, 5 2 , and C e describing the pertur- 
bations at the e stage entirely in terms of the output values of the functions defined at the 
% stage. 

Let us now turn to the e-m transition. Again, from the joining of h(rj), h'(r)), and ^(77) 
one derives 



n 



B 1 e~ iyi s + B 2 e iyi s*\ , 



q ^2 , ' ,-, , ,—i<l\ 



C e + — — [fix e^ 1 (s + 3y 2 e"* d ) + B 2 (s* + 3y 2 e id )\ , 



3a 2 6a 2 t/2 



(72) 



(73) 



C 2 



2^ 



5a 2 



C e + — {Sx e-^ [e~ id (-5 + 2y 2 2 - Qty 2 ) + iy 2 ] 



+B 2 [e ld (-5 + 2y\ + 6iy 2 ) - iy 2 ] } , 
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(74) 



where d — y 2 — yi = ^(V2 — Vi), s = 2i + (y 2 — 2i)e~ ld . Everything at the m stage is 
known as soon as the coefficients Ci, C 2 , and C m are known. They are expressed in terms 
of the coefficients Bi, B 2 , and C e attributed to the e stage. Since these numbers, in turn, 
are known implicitly in terms of the coefficients attributed to the % stage, we have linked the 
very beginning with the very end. 

Our next step is to impose the requirement (65) and to choose the comoving coordinate 
system at the m stage. From Eqs. (72) and (69) one can find 

CiD = -^ai{s[3ui(l - iyi) + «i] - s*[3ui(l + iyi) + v x }) , (75) 



where 



D = 2y\ + ?-^.\s*{\ + i Vl ) - s(l - i Vl )] 

= 2y\ + (2 - 7 ) [2 - (2 + ym ) cos d-(y 2 - 2 Vl ) sin d] . (76) 

Equation (75) says how to choose Cj at the i stage in order to match right to the comoving 
coordinate system at the m stage. We can now put Eq. (75) into Eq. (70) and solve Eqs. (70) 
and (71): 

bu ai 

B * ^ = -^-{^1 + <yi)«i - [(2 - 7 )a - 2yi]^i) • (78) 
bu ol\ 

Substituting these formulae and Eq. (75) into Eqs. (73) and (74), we reach our goal — the 
finding of Ci, and C 2 in the comoving coordinates: 

d = - ai {3u iyi (sin d + y x cos d) - uj[(2 - 7)(cos d - 1) - y x sin d]} , (79) 
6a 2 a\U 

C 2 = - Ql {3«iyi[(10 - 3y 2 - 10yiy 2 ) sin d + (-10y 2 + 10j/i - 3yiy 2 ) cos d] 
- Vl {-2(2 - 7 ) (5 + y 2 2 ) - fe/i(10 - 3t/ 2 2 ) - 10y 2 (2 - 7)] sin d 
+ [(2 - 7 )(10 - 3y 2 2 ) + Wy iy2 ] cos d}} . (80) 
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So far, no approximations have been made. We will start making them now. 

The numerical value of the denominator D, and hence the absolute values of C\ and C 2 , 
depend critically on how close the exiting value 7(771) = 71 is to 2. We are interested in 
wavelengths that were longer than the Hubble radius at 77 = r\\. Their wave numbers satisfy 
the requirement y± <C 1. On the other hand, yilvx = 02/^1 1, and d = y 2 — Hi ~ Vi- 
The wavelenghts longer than the present day A c = 2ira{r]R)/n c ~ 220 Mpc correspond to 
small z/2, and n < n c . These are the wavelenghts of the major interest for the discussion of 
the large-angular-scale anisotropy in CMBR. In the approximation of small 1/2, two leading 
terms in D are 



The second term is much larger than the first one (and, hence, the expected C±, and C2 are 
hopelessly small) for all 



unless the exiting value 7x is so close to 2 that the second term can be neglected. In order 
to deal with the most favorable situation and not proliferate complications, we will assume 
that 71 = 2 and (7 / /7)( r ?i) = 0- We know, see Sec. IV, that the transition from the very end 
of the % stage (after a thin "sandwich" interval) to the very beginning of the e stage can be 
made arbitrarily smooth (at least, in theory). So, we will be using D ks 2yf <C I. 

We should now take into account the fact that v 1 U\ for all wavelengths of our interest. 
Indeed, we are interested in modes that have interacted with the potential barrier in Eq. (38) 
and have been amplified at the % stage. For the scale factors (6), their wave numbers satisfy 
the requirement {nr]i) 2 /3(/3 + 1) which translates into the condition y\ -C /3/3(/3 + 1), or 
simply ?/j C 1. We will derive the approximate formulae valid in the leading order by the 
parameter y\. 

As we know, there are two independent solutions in the under-barrier region. Which of 
them dominates is the matter of choice of the initial conditions at i] = r\ (choice of constants 



D ss Ti2/i + 
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Al, and A 2 in Eq. (40)). For classical solutions, one can choose the initial data in such a 
way that there will be no amplification at all, or there will be even attenuation instead of 
amplification. However, a "typical" choice of initial data at r\ = r) , which amounts to the 
averaging over the initial phase (or a rigorous quantum- mechanical treatment), always leads 
to the dominant solution ji ~ 0^/7, and to amplification [19]. This is the choice that we 
imply here and will justify later, Eq. (102) in Sec. VIII. Since v ~ (fx/ay/j)' and ji ~ 
the quantity v\ is relatively small. Concretely, for solutions (40), one has 

M * 2 ^r 4 (^ + |) (nr?iri (81) 

and 

T) A 

W - H(Vi) « - I^r— ; 3, (2/3 + S)- 1 ^)^ 2 

2^+2 r(/3 + §) 

so that 

The finite jump of V\ while going through the "sandwich", see Sec. Ill, does not change this 
conclusion. Thus, we can neglect all the terms containing v\ in Eqs. (77)- (80). 
In the leading order, one has 

B 1 « -B 2 « v^fa) = 5 
fix + fi 2 « Byf . 

These formulae ensure the standing wave pattern for all wavelengths at the e stage and 
the power spectrum modulation at the m stage (compare with Sec. V). The leading order 
expressions for C±, C 2 are as follows: 

1 1 

Ci« \/2/x(r/i)sin d (82) 

2a 2 yi 



3 1 

C 2 ~ -7 v^/x^i) 

5a 2 yi 



'10 — 3y|) sin d — 10y 2 cos d 
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(83) 



We will give a brief analysis of Eqs. (82) and (83). The growing and decaying components 
of hi(rj) and Se/e , see Eqs. (66) and (67), are of the same order of magnitude at 77 = r\ 2 for 
all wavelengths. The coefficients C±, C 2 are smooth for long waves, d ~ y 2 <C 1, n <C n c : 



and are oscillating for shorter waves, n > n c (Sakharov' s oscillations). At a series of frequen- 
cies, the factor sin d is zero, and the growing component totally vanishes: no gravitational 
field perturbations, no time derivatives of the perturbations, no energy density perturba- 
tions. These modes were highly excited, like others, at the end of the % stage, but were 
stripped off of their energy by the very late times of their evolution. In terms of quantum 
mechanics, one can say that these modes have been desqueezed, sent back to the vacuum 
state [23]. The position of zeros is determined by ^(V2 — Vi) = ^h, k = 1,2,3 ... or, 
approximately, by ^- = nk. If one defines the distance x traveled by sound waves between 
the barriers at 77 = 771 and 77 = r\ 2 by x = a ^(V2 — Vi), the zeros arise when x is covered 
by an integer number of half-waves, x = The first zero in the spectrum of the growing 
component arises at k — 1 which corresponds to the present day scale of the order of 70 Mpc. 

The numerical values of C\ and C 2 as functions of n are controlled by the n-dependent 
function //(771). For simple scale factors (6) and solutions (40), /i(r/i) is given by Eq. (81) 
where the value of A 1 is determined by quantum mechanics, as will be discussed in Sec. VIII. 
The function /i(r)i) is exactly the same as the one used for gravitational wave calculations [17]. 
This allows us to make certain comparisons of density perturbations with gravitational 
waves. For instance, the growing component of gravitational waves taken at the beginning 
of the m stage, h gw (rj 2 ), has the following amplitude in the low frequency limit, n n c : 



This number is 3v / 3~7r times larger than Ci, Eq. (84), which gives the amplitude h for density 
perturbations in the same limit. We can also compare the growing component of hi with 
the growing component of gravitational waves. Let us take 77 = r] R and n n H where 




(84) 





vim) 



a/2 ai 
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riH = 47r/ (rjR — r) m ) corresponds to the wavelength equal to the Hubble radius at r\ = rjR. 
One can derive 



As we see, gravitational waves and density perturbations "enter" the (time-dependent) Hub- 
ble radius with approximately equal amplitudes, regardless of the numerical values of pa- 
rameters l , P describing the % stage. 

According to our definitions, see Sec. VIII, the Fourier component of the quantized 
field includes the factor l p i/y/2n in addition to h(r/) or hi{rf). We can give an estimate 
for the "characteristic" amplitude h(n) ~ nl v iC\ of the h field, which is a substitute for a 
more rigorously defined expectation value of the dispersion (square root of the variance) of 
the field. Combining Eqs. (84) and (81) we can find in the low frequency limit n <C n c : 
h(n) ~ (lpi/l )n l3+2 , i.e., exactly the same behavior as for gravitational waves. The growing 
components of ^ (77) and 5e/e (r)) contain the additional factor n 2 {rj — r] m ) 2 which gives two 
extra powers of n in their spectra. There is nothing spectacular about the De Sitter case 
(3 = —2. The derivative of the Hubble parameter can be arbitrarily close to zero at the 
time when the wavelength of our interest leaves the Hubble radius at the % stage. The 
perturbation will have a finite, not infinite, amplitude today. 




36 



VII. DENSITY AND ROTATIONAL PERTURBATIONS IN THE 
HIGH-FREQUENCY LIMIT 



The normalization of the perturbations is determined by quantum mechanics. We intend 
to amplify the zero-point quantum fluctuations of the primeval matter which is, in our case, 
the scalar field (1). Before the amplification, the frequencies of the fluctuations were much 
higher than the frequency of the gravitational pump field. To the modes of our interest 
the surrounding space-time seemed at the beginning almost flat. As a preparation for 
quantization, we will first consider density and rotational perturbations in matter placed in 
the Minkowski space-time, that is when gravity is totally neglected (a(rj) = 1 in Eq. (2)). 

The deformation of an elastic medium is usually described [24] with the help of a dis- 
placement three-vector Wj(t,x) which can be written as 

Ui = £(t)Q,i + 0(t)Qi ■ (85) 
The scalar function Q is defined by Eq. (8), the vector function Qi is defined by the equations 

Q l , k ' k + n 2 Q t = 0, Q\ l = . (86) 

The deformation tensor is 

11 1 

u ik = -j{ u i,k + Uk,i) — -ji{Q,i,k + Q,k,i) + ^(Qi,k + Qk,i) 

and its trace is u — u l ti = —^n 2 Q. The stress tensor can be written in the general form 

= s(Q :ijk + Q :k:i ) + pc 2 6(Q ijk + Q kji ) + n 2 (2s - pc 2 £)Q5 ik (87) 

where p is density, q is the longitudinal velocity of sound, c t is the transverse (torsional) 
velocity of sound, and s is arbitrary function of time. The equations of motion 

d 2 Ui da\ 
^ dt 2 dx k 

reduce to the oscillatory equations for elastic waves. In terms of 77-time, they can be written 

as 
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£" + n 2 ^ = 0, 6" + n 2 ^6 = . (88) 
c c 

We shall now relate the theory of elasticity with the theory of cosmological perturbations. 
We shall consider the high-frequency limit of the perturbations, that is when the scale factor 
a (r/) is almost constant and its variability can be neglected in comparison with frequencies 
of the waves. As we know, the perturbed components of the energy-momentum tensor for 
density and rotational perturbations have the general form (see Eq. (10) and Ref. [2]): 

T =-^iQ T° = -T* = ^'Q, + ^0'Q*, 
a 2 a 2 a 2 

T '= + - hH Qt * + «*•<) + ■ (89) 

The differential conservation laws T@ „ = can be reduced in the high-frequency limit to 
the following equations: the a = component gives the equation — e[ + £'n 2 = 0, which 
integrates to 

n\ = e 1 , (90) 

the a — i components give 

£"+Pi = 0, 9" + X = . (91) 

Equations (90) and (91) can, of course, be obtained from the perturbed Einstein equations 
as well. 

The stress tensor is connected with the perturbed components by = —pc 2 Ti k 
(a(rj) = 1). From the comparison of Eqs. (87) and (89) one finds 

f.=" 2 | { - x=n2 i e - s = -" c2 |j- (92) 

With these expressions for pi, x, Eqs. (91) coincide with the wave equations (88). In 
cosmology and theory of elasticity, we are dealing essentially with the same physics. 

The quantization of density and rotational perturbations should be based on Eqs. (88). 
Rotational perturbations have been considered elsewhere [2]. The quantum-mechanically 
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generated rotational perturbations can contribute to the CMBR anisotropy and may be 
important for the smaller scale astrophysics. (One should be aware, though, that there 
exists also an alternative view on the subject according to which rotational perturbations 
are "irrelevant for cosmology" [25].) We will concentrate here on density perturbations. 
For simple models of matter, such as perfect fluids and scalar fields, the functions x an d 
pi vanish. One is left with the single variable £ and isotropic stresses. The quantization of 
these oscillations in an elastic material placed in the Minkowski space-time would lead to 
the notion of phonons. 

There is no wonder that in the case of scalar field matter the role of £ is played by <pi. 
If one writes £(77) = £ e~ mri and tpi(r)) = ip w e~ mri , Eq. (24) gives 

in£ = tp' Q tpva . (93) 

In the high-frequency limit (large n), the first term in Eqs. (23) and (25) dominates. With 
the help of Eq. (93) one derives e\ — p\ — n 2 £. In other words, for the high-frequency 
scalar field perturbations, the velocity of sound is almost equal to the velocity of light. The 
perturbations behave as massless scalar particles. 

It is the scalar field oscillations that should be normalized by ascribing a "half of the 
quantum" to each mode. Due to the Einstein equations the scalar field perturbations are 
accompanied by the gravitational field perturbations. In the high-frequency limit, the second 
term in Eq. (43) can be neglected, the normalization of ipi transfers to the gravitational field 
variable ji and ultimately to h and hi. This is how we will know the initial amplitude for 
density perturbations. 
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VIII. QUANTIZATION OF DENSITY PERTURBATIONS 



In the limit of a free massless scalar field placed in the Minkowski space-time, we have 
for each mode: 

S Vk = Vl Q + V IQ* = Vl (t) e iky + <pl(t) e~ iky . 
The total field can be written as 

I r°° 1 r i n 

M*. y) = # k7 = [* e-« e * + 4 e« e"<*] . (94) 

The normalization constant C is to be found from the requirement 

/OO j ("OO 

ed 3 y|0> = ^ / d 3 kw fe (0|c k c k + 4c k |0) , 
-OO Z J —OO 

where e is the energy density of the field. Since in our case 

1 



6= 2 



we derive C = cy/k. Obviously, the normalization coefficient C includes the Planck constant 
h but does not include the gravitational constant. 

To write the field operator (94) in the curved space-time (2) one should make the following 
replacements: 

y = a(?7)x, k=— !— n, w k = , c k = [a(^)] 3/2 c n . 
a{rj) a(rj) 

The field operator takes on the form 

1 1 POO 1 

6<p(r,, x) = -tt^--. / rf 3 n-=[c n e-*™» e* nx + <£ e *™> e— x ] . (95) 
a{rj) J-oo y/2n 

This expression is only valid in the high-frequency limit, when the field can be regarded as 
free and its non-adiabatic interaction with gravity can be neglected. When the interaction 
becomes important, the time dependence of the field ceases to be so simple. The operator 
c n e~ mri , should be replaced by c n (rj), and the evolution of c n (i]), c^irf), should be found from 
the Heisenberg equations of motion. 
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The matter perturbations are accompanied by the gravitational field perturbations. Due 
to the Einstein equations they are linked together and form, in a sense, a united entity. As 
we have seen in Sec. Ill, the entire dynamical problem at the i stage reduces to a single wave 
equation (38) for a single variable /i(r)). All perturbations can be found from a given solution 
to this equation. It follows from Eq. (43) that \^2~Kaipi(r]) in the high-frequency limit. 

Since the positive frequency scalar field n-mode solution is <pi(rj) ~ (l/a(r]))\^chc n e~ m ' n 
this leads to /j,(r}) ~ Ay/ir l p ic n e~ mri . Note that the normalization coefficient includes the 
gravitational constant and is proportional to the Planck length l p i = (Gh/c 3 ) 1 ^ 2 . 

Having in mind the basic wave equation (38), we can now introduce the "fundamen- 
tal" scalar field $(77, x) which describes the whole quantum system interacting with the 
gravitational pump field: 

roo 1 



1 r°° 



Pn(»7)e il,x + ctfo)e 



(96) 



(27T) 3 / 2 J-oo v 7 ^ 

The annihilation and creation operators c n (r]), c^rj) are governed by the Heisenberg equa- 
tions of motion 

^M = - iM ,),tf], *M = - ¥M . (97) 

The interaction Hamiltonian H is given by 

H = nclcn + nc^ n c_ n + 2a(r])cl 1 cl n + 2a* (?])c n c_ n , (98) 

where the coupling function a{rf) is <j(i]) = | . Equation (98) demonstrates explicitly 
the underlying parametric interaction of the field oscillators with the pump field. This is 
simply a generalization of a theory previously developed for gravitational waves (for a review, 
see Ref. 19). 

The common way of solving Eqs. (97) and (98) is to write the operators in the form 
cM = u n ( V )c n (0) + v^cl^O) , (4(7?) = <(77)4(0) + <(l)c-„(0) , (99) 

where c n (0), 4,(0) are the initial values of the operators taken at some i] = t/q long before the 
interaction became effective, and [c n (0), c^(0)] = <5 3 (n — m) The classical complex functions 
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Un{v), v n(v) (do not mix up with the functions u(rj), v{rj) introduced in Sec. Ill) obey the 
condition |w n | 2 — \v n \ 2 = 1 and satisfy the equations 

iu' n = nu n + 2av* a , iv' n = nv n + 2au* n (100) 

with the initial data u n (0) = 1, v n (0) = 0. If one introduces /i n (r)) = u n (r)) + v*(r)), it 
follows from Eqs. (100) that the function fi n (v) should satisfy precisely the Eq. (38). The 
initial conditions for [i n (v) m the high-frequency limit \nr]\ — > oo are fJ, n (v) ~^ e~ m ^~ Vo \ 
Vniv) -ine~ in ^ a) . 

For each mode n there exists the vacuum state |0 n ) defined by the condition c n (0)|0 n ) = 0. 
As a result of the Schrodinger evolution with the Hamiltonian (98), the initial vacuum state 
|0 n ,-n) = |0 n )|0_ n ) transforms into a multiparticle two-mode squeezed vacuum state (see 
Ref. 19 and references cited therein). In other words, the perturbations (waves) are generated 
in correlated pairs. The statistical properties of the field are determined by c n (r]), cj,(r/). By 
using Eq. (99) one can rewrite Eq. (96) in the form 

1 Z" 00 1 



j roc 

$(^,x)=4^— — J <Pn- 



c n (0)/x n (r / )e- x + 4(0)/x;(r / )e- 



(101) 



(2tt) 3 /2 J_ x ^ 

where the functions /i n (^), Univ) should be taken with the appropriate initial conditions 
discussed above. For simple solutions (40), the initial conditions translate into the require- 
ments 

A 1 = 1 .fl f>(™*+¥) A 2 =iA ie - in P . (102) 

cos pir V 2 

The quantized gravitational field perturbations are expressible entirely in terms of the 
$(?7,x) field (96). There are many components of hij but there is only one sort of creation 

1 2 

and annihilation operators. Let us introduce new notations h(i]) = h(r)), h{f]) = hi(rj) and 

1 2 1 2*3 

the polarization tensors = 5^-, = —nirij/n , and PijP = — 1. The field operator 
hij(r),yi) can be written as 

I roc 1 2 



I f°° ! s 

h i3 = 4V5FZ„^ d^n-^ g P y (n) 



c„(0)/ in e inx + cj 1 (0)/ ln e- 



(103) 
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where the classical complex functions h n (v) should be derived from n n {rj) through the equa- 
tions and initial conditions already discussed. A similar expression can be written for the 
operator of the energy density perturbation 5e/e . Equation (103) is the starting point for 
the calculation of the expected angular anisotropy in CMBR. 
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IX. VARIATIONS OF THE CMBR TEMPERATURE CAUSED BY DENSITY 
PERTURBATIONS OF QUANTUM-MECHANICAL ORIGIN 



The photons of CMBR are emitted at r] — rj E and are received by us at r] = rj R . A 
particular direction of observations is characterized by the unit vector e k = (sin 9 cos <j>, 
sin 9 sin 0, cos 9). In absence of perturbations, temperature of CMBR seen in all directions 
would be the same, T. Gravitational field associated with the density perturbations at 
the m stage causes a variation of the temperature with respect to the unperturbed value 
T [26]: 

»ifdhi 



ST k 1 r^fdhij t A, 



(104) 



where w\ = r) R — tje, and dhij/drj is taken in the comoving synchronous coordinate system 
along the path x k = e k w, rj — t)r — w. In case of small perturbations, which we are actually 
dealing with, the emission time r] E can be regarded as being one and the same for all 
directions. Since the scale factor satisfies the approximate relationship a(rjE) /o,(tIr) ~ 10~ 3 , 
and the rj time can be chosen in such a way that r] R — r\ m = 1, the quantity W\ is close to 1. 
We will use w\ — l. The wavelength equal to the present day Hubble radius Ih corresponds 

to TlH = 47T. 

For the quantized perturbations, the temperature variation 5T/T, Eq. (104), becomes 
a quantum-mechanical operator. Since it is only the longitudinal-longitudinal part of hij 
that participates in producing 5T/T, we can write 



ST 

T 1 



e k ) = -2y/^l 



1 



P i 



(2tt) 3 / 2 Jo 



dw 



-oo n 2 



x 



c n (0)/ n ( % - w)e^ e w + ci(0)r n ( VR - w)e- in ^ 



(105) 



where 



fn(VR -w) = 



1 dhtirj) 



2n drj 



The individual observed distributions of the CMBR temperature over the sky should 
be compared with theoretical predictions based on the quantum-mechanical expectation 
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values. The mean value of 5T/T(e k ) is obviously zero, (0\5T/T (e fc )|0) = 0. The variance 
(0\5T/T (e k )5T/T(e k )\0) is not zero but does not depend on the point and direction of 
observations. To study the angular distribution of the temperature variations one should 
construct the angular correlation function K for two different directions e k and e k , 

e-i^ki = cos 5: 

K = (0\^(e k )^(e k 2 )\0) . 
By manipulating with the product of two expressions (105), one can derive 

I rl rl roc rn r2n; (ri-p^X^ (n-p^X^ 

K = ^l 2 pl -^—l dw dw n 2 \f n \ 2 dn sm9d9 l^iJ ^e 2 j cos(wfcC fc )d0> 

[ZiTT ) J JO JO JO JO Tt 71/ 

(106) 



where ( k = e\w — e k w. A lengthy calculation of the integrals over the angular variables 0, 
9 gives the following result: 

r*n6d0 1^ cos(n fc C fc )# 
Jo Jo n 2 n z 

= An^cos 2 5(nCr l/2 J 1/2 (nC) + (1-5 cos 2 5)«)- 3 / 2 J 3/2 «) 

+2 cos 5( 1 - cos 2 5) (nw) (nw) (n()~ 5/2 J 5/2 «) +4(3 cos 2 5-1) «)~ 5/2 J 5/2 (n() 
+8 cos 5 (cos 2 5 -l)(nw)(nw)(n()- 7 / 2 J 7/2 (n() + (cos 2 5-l) 2 (nw) 2 (nw) 2 (n()- 9/2 J 9/2 (n()} , 

(107) 

where £ = (w 2 — 2wwcos5 + w 2 ) 1 ^ 2 . For further calculations one may use the following 
formula (valid for half-integer v): 

(n() U<) = v^J> + *0 {nw)v {nnj)u dzV _ 1/2 P k+v - 1/2 (z) , (108) 

where z = cos 5 and Pi(z) are the Legendre polynomials. (I derived and used this formula 
in course of studying the gravitational wave [17] and rotational [2] perturbations, but I 
believe that this formula may exist somewhere in the previously published literature.) With 
the help of Eq. (108), one can rearrange the correlation function K to the following final 
expression: 
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K = l 2 pl Y^K l P l {cos5) , 



(109) 



1=0 



where 



K t = (21 + 1) / \I ln \ 2 dn 
Jo 



= f 

JO 



fn(VR ~ x/ri) 



X 



1/2 



1(1-1) 



j h+l (x)+h l+l (x)y x , (no) 



and x = nw. The derived formula for K is general and can be used with arbitrary function 
/„. In practice, the limits of integration over n are determined by the frequency interval 
within which the perturbations were really generated. 

As we see, the decomposition of K consists of all multipoles including the monopole, 
I = 0, and dipole, I = 1, terms. To carry out the calculations up to a concrete number, we 
will consider scale factors (6) and solutions (40). 

As we know, the growing and decaying components of hi(rj) are of the same order of 
magnitude at i] = i] 2 . However, the decaying component is decreasing since then and can 
be neglected in the calculation of 5T/T. For the coefficient C\ responsible for the growing 
solution, see Eq. (82), we have 

1 



\Ci\ 



211 



2 2/3 2-2 

n p n c sin 



n 



where 



(111) 



7T 



2^cospnT(p+ -) 
2 



1 for/3 = -2 



The expression for Ki takes on the form 

K x = -±-\m\ 2 n 2 c (2l + 1) rn^W-(i,„) 2 dn, 
100/q jo n c 

where 



Hn = 



n — x 



o x 



1/2 



1(1-1) 



r 



Ji+i(x) + -J3+i(x) \dx 



(112) 



(113) 
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We will now estimate the contribution of long waves, n < 1, to the lower order multipoles 
K t . The integrals % Xn should be calculated separately for I — 0, I — 1, and I > 2, because 
of the factor 1(1 — 1) in Eq. (113). For I > 2, the term with this factor dominates. The 
approximate expression for (i/„) 2 , I > 2, takes on the form 



For / = and I = 1, the term with the factor 1(1 — 1) does not contribute to (ii n ) 2 . The 
result still have the form of Eq. (114) but n 21 should be replaced by j^n 6 for I — 1, and by 
^n 4 for I — 0. These results are in full agreement with Ref. 13: in the limit of long waves, 
the monopole and dipole contributions of an individual wave are suppressed; the monopole 
component is of the same order of magnitude as the quadrupole component, while the dipole 
component is further suppressed by an extra power of n. 

We should now use (ii n ) 2 for the calculation of K t . This is where the spectrum of the 
perturbations comes into play. We can write for I > 2: 



and the appropriate replacements discussed above should be made for I — 1, I — 0. Since we 
are working in the limit of small n, the integration over n cannot be extended to the values 
n > 1. However, typically, short waves contribute little to the lower index multipoles. The 
values of n up to n ~ 21 are more important, but for the purposes of simple evaluation we 
restrict the integration by n = 1. 

The suppression of the monopole contributions of individual waves saves us from a big 
trouble. If Eq. (115) were true for / = 0, the monopole term K would be power-law divergent 
in the limit of n — > for all (3 < —2. In order not to be in conflict with the finite observed 
2.7 K temperature, we would need to resort to the fine tuned minimally sufficient duration of 



{knf » A' 




(114) 



where 





(115) 



47 



the i stage, in which case the long waves with this spectrum are simply not being generated. 
From the correct Eq. (115) follows that the danger of divergence for K and K 2 arises only 
in models with (3 < —4. (The quadrupole anisotropy produced by gravitational waves does 
also diverge in these models [17].) Thus, the interval — 2 > (3 > — 4 is potentially allowed. 

We will now introduce the notations l v i\[K~§ = M, l v i^/K[ = D, l p i\fK~2 = Q, and will 
compare M, D, and Q. For the quadrupole Q, one can find from Eq. (115) 

Q "^^W7l ' (116) 

The monopole M and dipole D are related with Q by 



M _ y/E D 13 1(3 + 4 
Q~~1T' Q ~ V 20 V P + 5 ' 

We do not have an observational access to the unperturbed temperature T, but whatever 
is the measured Q, we can expect that a correction of about the same magnitude as Q is 
included in the measured T. The same is true for the dipole component D (there is little 
doubt, however, that the overwhelming part of the measured dipole anisotropy is accounted 
for by our peculiar motion). 

We will now compare the contributions of density perturbations and gravitational waves 
to the components K\ of the correlation function K in the long wave limit, n < 1. We should 
compare Eq. (115) with the analogous expression for gravitational waves [17]. The ratio of 

9 d 

the gravity wave contribution Ki to the density contribution Ki has the form 

Ki (/ + l)(/ + 2)(2/ + l) 2 



ki 21(1-1) 

We see that the ratio is independent of the parameters l a , (3 describing the % stage. For the 

9 d 

quadrupoles Q and Q, we have in the long wave limit 

9 

Q i— 
Q 

that is a somewhat larger contribution of gravitational waves. It is necessary to take into 
account also the shorter waves in order to get a more accurate estimate. 
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In conclusion, there is no dimensionless ratios that could be adjusted in such a way 
that the contribution of density perturbations to the quadrupole anisotropy would be much 
larger than the contribution of gravitational waves. These contributions are of the same 
order of magnitude while numerical coefficients are somewhat in favour of gravitational 
waves. At the same time, the very generation of density perturbations (and rotational 
perturbations) is more problematic than the generation of gravitational waves. On these 
grounds one can conclude that if the observed large-angular-scale anisotropy of CMBR is 
caused by cosmological perturbations of quantum- mechanical origin (what else?), they are, 
most likely, gravitational waves. 
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